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Abstract The hadronic corrections to the muon anomalous 
magnetic moment a u , due to the gauge-invariant set of di¬ 
agrams with dynamical quark loop light-by-light scattering 
insertions, are calculated in the framework of the nonlocal 
chiral quark model. These results complete calculations of 
all hadronic light-by-light scattering contributions to cq, in 
the leading order in the 1 / N c expansion. The result for the 
quark loop contribution is a |J LbL ' Loo P _ (11.0 ±0.9) • 10 10 , 
and the total result is a |) Lb k N *Q M _ g ±1.2)-10~ 10 . 

1 Introduction 

Experimental and theoretical research on lepton anomalous 
magnetic moments has a long and prominent history * 1 . The 
most recent and precise measurements of the muon anoma¬ 
lous magnetic moment tq, were published in 2006 by the 
E821 collaboration at the Brookhaven National Laboratory 
[5]. The combined result, based on nearly equal samples of 
positive and negative muons, is 

= H6 592 08.0 (6.3) x 10“ 10 [0.54 ppm], (1) 

Later on, this value was corrected [6, 7] for a small shift 
in the ratio of the magnetic moments of the muon and the 
proton as 

a BNL,CODATA = jjg 592 Qq j (53) x JO -10 . ( 2 ) 

This exiting result is still limited by the statistical errors, 
and proposals to measure cq, with a fourfold improvement in 
accuracy were suggested at Fermilab (USA) [8] and J-PARC 
(Japan) [9]. These plans are very important in view of a very 

a e-mail: dorokhov@theor.jinr.ru 
b e-mail: aradzh@icc.ru 
c e-mail: zhevlakov@phys.tsu.ru 

1 For comprehensive reviews see [ 1—4], 


accurate prediction of cq, within the standard model (SM). 
The dominant contribution in the SM comes from QED 

= 116 584 71.8951(80) x 1(T 10 [10]. (3) 

Other contributions are due to the electroweak corrections 

[ 11 , 12 ] 

a™ = 15.36(0.1) x 10~ 10 [12], (4) 

the hadron vacuum polarization (HVP) contributions in the 
leading, next-to-leading and next-next-to-leading order [13- 


15], 




HVP.LO 

= 

694.91(3.72)(2.10) x 1(T 10 

[13], 

(5) 

HVP,NLO 
a p — 

—9.84(0.06)(0.04) x 1(T 10 

[13], 

(6) 

HVP.NNLO 

= 

1.24(0.01) x 10~ 10 

[14], 

(7) 


and the hadronic light-by-light (HLbL) scattering contribu¬ 
tion (as it is estimated in [16]) 

^HLbL _ 10.5(2.6) • 10~ 10 . (8) 


As a result, the total value for the SM contribution, if we 
take (8) for HLbL, is 

a® M = 116 591 84.1 (5.0) x 10~ 10 . (9) 

From the comparison of (2) with (9) it follows that there 
is a 3.11 standard deviation between theory and experiment. 
This might be an evidence for the existence of new inter¬ 
actions and stringently constrains the parametric space of 
hypothetical interactions extending the SM 2 . 

From above it is clear, that the main source of theoret¬ 
ical uncertainties comes from the hadronic contributions. 

2 In this regard we would like to mention the work [17], where a public 
code for computing new physics contributions to a^ applicable to any 
particle physics models is developed. 
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The HVP contribution flJJ VPLO , using analyticity and uni- 
tarity, can be expressed as a convolution integral over the 
invariant mass of a known kinematical factor and the total 
e + e~ —> y* —> hadrons cross-section [18]. Then the corre¬ 
sponding error in «| 1 IVI> L<> essentially depends on the accu¬ 
racy in the measurement of the cross-section [13, 15]. In 
near future it is expected, that new and precise measure¬ 
ments from CMD3 and SND at VEPP-2000 in Novosibirsk, 
BES III in Beijing and KLOE-2 at DAFNE in Frascati will 

allow to significantly increase the accuracy of the predic- 
.- f HVP,LO 
tions for a^ 

On the other hand, the HLbL contribution aJ) LbL cannot 
be calculated from first principles or (unlike to HVP) di¬ 
rectly extracted from phenomenological considerations. In¬ 
stead, it has to be evaluated using various QCD inspired 
hadronic models that correctly reproduce basic low- and high- 
energy properties of the strong interaction. Nevertheless, as 
will be discussed below, it is important for model calcula¬ 
tions, that phenomenological information and well estab¬ 
lished theoretical principles should significantly reduce the 
number of model assumptions and the allowable space of 
model parameters. 

Different approaches to the calculation of the contribu¬ 
tions from the HLbL scattering process to aJJ LbL have been 
suggested. These approaches can be classified into several 
types. The first one consists of various extended versions of 
the vector meson dominance model (VMD) supplemented 
by the ideas of the chiral effective theory, such as the hid¬ 
den local symmetry model (HLS) [20], the lowest meson 
dominance (LMD) [21-23], and the (resonance) chiral per¬ 
turbative theory ((R)XpT) [24—26]. The second type of ap¬ 
proaches is based on the consideration of effective models 
of QCD that use the dynamical quarks as effective degrees 
of freedom. The rest include different versions of the (ex¬ 
tended) Nambu-Jona-Lasinio model (E)NJL [27], the con¬ 
stituent quark models with local interaction (CQM) [28- 
32], the models based on nonperturbative quark-gluon dy¬ 
namics, like the nonlocal chiral quark model (NjQM) [33- 
39], the Dyson-Schwinger model [40] (DS), or the holo¬ 
graphic models (HM) [41, 42]. More recently, there have 
been attempts to estimate within the dispersive ap¬ 

proach (DA) [43, 44] and the so-called rational approxima¬ 
tion (RA) approach [45]. 

The aim of this work is to complete calculations of the 
leading in 1 /N c HLbL contributions within the N^QM started 
in [38, 39] and compare the result with (8). Namely, in pre¬ 
vious works we made detailed calculations of hadronic con¬ 
tributions due to the exchange diagrams in the channels of 
light pseudoscalar and scalar mesons. In the present work, 
the detailed calculation of the light quark loop contribution 
is given 1 . 


Preliminary results of this work were announced in [3], 


2 Light-by-light contribution to in the general case 

We start from some general consideration of the connection 
between the muon AMM and the light-by-light (LbL) scat¬ 
tering polarization tensor. The muon AMM for the LbL con¬ 
tribution can be extracted by using the projection [46] 

a ]l hL = + Y}{p + m ^U pa (p,p)) , 

( 10 ) 


where 


tt ( ' \ - 6 l d4qi l d4qi _ \ _ 

pa P 11 J (2n) 4 * J (2n) 4 q\(q x + q 2 ) 2 (q\ + k) 2 


,...u P'-qi+m M p + q i+m M , 

Y (p'-q 2 ) 2 -ml Y (p + qi ) 2 -ml Y 


X 

( 11 ) 


x -^^^via(q2,-(qi+q2),k + qu-k), 


where is the muon mass, = (p 1 — p) fll and it is neces¬ 
sary to make the static limit k^ —> 0 after differentiation. Let 
us introduce the notation 


-g^- u nvia(q 2 ,-{qi +qi),k+qu-k) = 

n pnvxa(q2,-(qi + q2),qi) + 0(k) (12) 

for the derivative of the four-rank polarization tensor 4 , and 
rewrite Eqs. (10) and (11) in the form (q 2 = q x + q 2 ) 


LbL e 6 f A1 f A2 ,, 

p 48 J (27t) 4 J (2n) 4 

J Jp/ivAg(g2,-g3,gi)T p/n ' ;l0 ' (q\,q 2 ,p) 

q\qlq\{{p+q\) 2 - ml)((p - q 2) 2 - m 2 ) ’ 


where the tensor '\'Pt- lvla is the Dirac trace 


T p ^ vXa ( quq 2 , P )=Tr(^p + mil )[f,Y a ]{p + m M ) x 
x y p (p-q 2 + m^)f (p + q x +w M )y A 

Taking the Dirac trace, the tensor r j'PP v2 - a becomes a 
polynomial in the momenta p, q\, q 2 . After that, it is conve¬ 
nient to convert all momenta into the Euclidean space, and 
we will use the capital letters P, Q \, Qp for the correspond¬ 
ing counterparts of the Minkowskian vectors p, q\, q 2 , e.g. 

4 First note, the tensor Ylp V x a can be of any nature (QED, hadronic, 

etc.) Another note concerns the important result expressing the tensor 
Tlp^vXa in the explicitly gauge-invariant form that was obtained in [47, 

48], 
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P 2 = -p 2 = -m 2 , Q\ = ~q\, Q\ = -q\. Then Eq. (13) be- 
comes 

LbL _ e 6 f dpQl I" d%Q 2 1 TP^ vA(T n pMvAg 
P 48m p J (2n) 4 J (2n) 4 Q\Q\Q 2 D x D 2 

D x = {P + Q x f + ml=2(P-Q l ) + Q\,x (14) 

!h = (P~ Qi? + ml = -2(P • Q 2 ) + Q 2 . 


Since the highest order of the power of the muon mo¬ 
mentum P in XPP vA<J is two 5 and TI ppv2ia is independent of 
P, the factors in the integrand of (14) can be rewritten as 


t p ^ a n ppvla 

D\P>2 


6 


= 2> a n„, 

a=l 


(15) 


with 

_ / 4m 2 

x=V^ 2 , Ri = Jl + -^(i=l,2), (19) 

z=^ (1 -^t )(1 -^ 2) - 

After averaging the LbL contribution can be represented in 
the form 

oo oo 

a LbL = JdQ x JdQ 2 p LbL (QuQ 2 ), (20) 

0 0 

with the density p LhL (Qi,Q 2 ) being defined as 


with the coefficients 


A 



A 2 


a 5 


1 

D^D 2 



A 6 = l, 


{P-Qi) 

D i 


(P^ 

P>2 


(16) 


where all P-dependence is included in the A a factors, while 
f[ a are P-independent. 

Then, one can average over the direction of the muon 
momentum P (as was suggested in [1] for the pion-exchange 
contribution) 

[ 4Qi f 4Q2 A a 

J (2n) 4 J {2nY Q\QlQl- 

oo oo 1 

^ J dQi f dQ 2 J dt fiQljk(A a )..., (17) 

00-1 

where the radial variables of integration Q\ = \Q\ | and Q 2 = 
\Q 2 \ and the angular variable t = (Qi ■ Q 2 )/ (|Gi| IQ 2 I) are 
introduced. The averaged A a factors are [1] 


<*>! = 


D, 


<A)2 = i \51 


P 2-1 

2 m 2 p 


m = { ^)= +{Ql . Ql) P^r 


D, 


8 m 2 p 


( a ) 4 = = -(< 2i - 22)- (1 ~ Rl)1 


{A) 5 = 


D 2 

DiD 2 / m 2 QiQ 2 x 


8m 2 p 


arc tan 


zx 

1 —zt 


(A) 6 = ( 1) = 1, 


(18) 


5 The possible combinations with momentum P are 


(p. QP 2 = (p. Q X )(D , - Q\)/2, (P- Q 2 ) 2 = -(P-Q 2 )(D 2 - Ql)/2, 

(P- Qi){P-Qi) = -Pi - Ql)(D 2 - 2l)/4, 

(P- Qi) = Pi “ Qi)/ 2, (P- fi>) = -(P2 - Ql)/2- 


p LbL (Gi,G 2 ) = ^ t [dt^^(A a )fl a . (21) 

27r " o=l_ A G 3 

Thus, the number of momentum integrations in the orig¬ 
inal expression for (10) is reduced from eight to three. The 
transformations leading from (10) to (20), are of general 
nature, independent of the theoretical (model) assumptions 
on the form of the polarization tensors fl„. In particular, 
this 3D-representation is common for all hadronic LbL con¬ 
tributions: the pseudoscalar meson exchange contributions 
[1, 38], the scalar meson exchange contributions [39], and 
the quark loop contributions discussed in the present work. 
The next problem to be elaborated is the calculation of the 
density p HLbL (Gi ■ G 2 ) in the framework of the model. 


3 Hadronic Light-by-light contribution to a p within 
N*QM 


Let us briefly review the basic facts about the N^QM 6 . The 
Lagrangian of the SU (3) nonlocal chiral quark model with 
SU (3) x SU (3) symmetry has the form 


% = q{x){id - m c )q{x) + — [A ?( x ) J s W 


H 

- JTabc 


/sM/sW/sM ~ 3Js(x)Jp S (x)Jp S ( x ) 


( 22 ) 


where q (x) are the quark fields, m c {m CyU = m c ^ ^ m c s ) is 
the diagonal matrix of the quark current masses, and G and 
H are the four- and six-quark coupling constants. The non¬ 
local structure of the model is introduced via the nonlocal 
quark currents 

Jm( x ) = J d 4 x x d 4 x 2 F(xi,x 2 )q(x — xi)r^q(x + x 2 ), (23) 


where M = S for the scalar and M = PS for the pseudoscalar 
channels, L" = A", r ps = iy 5 A", and P(xi,x 2 ) is the form 

6 More detailed information about the model is contained in our previ¬ 
ous works [37, 39], 
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factor with the nonlocality parameter A reflecting the non¬ 
local properties of the QCD vacuum. The SU( 2) version 
of the N^QM with SU (2) x SU (2) symmetry is obtained 
by setting H to zero and taking only scalar-isoscalar and 
pseudoscalar-isovector currents. 

Within the N^QM, the standard mechanism for spon¬ 
taneous breaking of chiral symmetry occurs, which is typi¬ 
cal for the Nambu-Jona-Lasinio type models with the chiral 
symmetric four-fermion interaction (local or nonlocal). Due 
to this interaction the massless quark becomes massive, and 
in the hadron spectrum the gap between the massless (in the 
chiral limit) Nambu-Goldstone pion and the massive scalar 
meson appears. This feature is common for the models used 
for the calculation of the hadronic contributions to the muon 
g — 2: the extended NJL model [27], the constituent chi¬ 
ral quark model [32], the Dyson-Schwinger model [40], the 
nonlocal chiral quark model [33-39]. In the nonlocal models 
the dynamically generated quark mass becomes momentum 
dependent and the inverse dynamical quark propagator takes 
the form 


9i(m) <?i(/i) q^(y) 




(c) (d) 

Fig. 1 The quark-photon vertex ( q ), the quark-two-photon vertex 

(?l, 92 ). the quark-three-photon vertex Fpyp(qi,q 2 ,q}), and the 
quark-four-photon vertex j 92 j 93 i 94 ) • 


S~ l (k) =k-m(k 2 ) (24) 

where m(k 2 ) = m c + moF(k 2 ,k 2 ) is the dynamical quark 
mass obtained by solving the Dyson-Schwinger equation. 
The significant feature of the nonlocal models [33-35] is 
that they correctly interpolate between the low-energy re¬ 
gion (and consistent with the low-energy theorems) and the 
high-energy region (where they are consistent with OPE). 
The basic fact is that the momentum dependent dynami¬ 
cal quark mass, that is the constituent quark mass m( 0) = 
m c + mo at low virtualities, becomes the current quark mass 
m c at large virtualities. This is in contrast to the local models, 
where the quark mass is the constituent one at any virtuality. 

For numerical estimates two versions of the form factor 
(in momentum space) are used: the Gaussian form factor 

E G (4,4)=exp(—4/A 2 ), (25) 


Then, apart from the kinetic term, the additional, nonlocal 
terms in the interaction of quarks with the gauge field are 
generated via substitution 

Jm{x) -f Jm(x) = jd 4 x\d 4 x 2 f{x\)f{x 2 )x 

x Q(x-xi,x)r^Q(x,x + x 2 ), (28) 

inducing the quark-antiquark-n-photon vertices. In order to 
obtain the explicit form of these vertices, it is necessary to 
fix the rules for the contour integral in the phase factor. The 
scheme, based on the rules that the derivative of the contour 
integral does not depend on the path shape 

y y 

Jdz v F v (z) = <5 (4) (x-y) Jdz v F v (z) = 0, 

X X 


and the Lorentzian form factor 

ft<44) = (7TW <26) 


The second version is used in order to test the stability of the 
results to the nonlocality shape. 

Next, it is necessary to introduce in the nonlocal chiral 
Lagrangian (22) the gauge-invariant interaction with an ex¬ 
ternal photon field A fl (z) . This can be made through the 
introduction of the path-ordered Schwinger phase factor for 
the quark field as 


q(y) Q{x,y) = ^expi^i£ dz tl A ll (z) 1 >q(y). 


(27) 


was suggested in [49] and applied to nonlocal models in 
[50] . For our purpose, we need to consider the quark-antiquark 
vertices with one-, two-, three- and four- photons (Fig. 1). 
The first two types of vertices was derived in [50], the ver¬ 
tex with three photons was obtained in [51], and the quark- 
4-photon vertex is given in the present work. Their explicit 
form and the definition for the finite-difference derivatives 
m <n> (k. k' ) are presented in the Appendix. The simplest quark- 
photon vertex has the usual local part as well as the non¬ 
local piece in terms of the first finite-difference derivative 

mP-\k,U) 

r^( qi ) = Y,+Ar^( qi ), 

4r ^ (1) (<7t) = - (k + k')^m {l) (k,k!) , 


(29) 

(30) 
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(a) (b) 


(c) (d) (e) 


Fig. 2 The box diagram and the diagrams with nonlocal multiphoton interaction vertices represent the gauge invariant set of diagrams contributing 
to the polarization tensor n fJv x p (<7i ,^2 i? 3,94)- The numbers in front of the diagrams are the combinatoric factors. 


while the quark-antiquark vertices with more than one pho¬ 
ton insertion are purely nonlocal (see Appendix). 

We have to remind, that in the models with the chiral 
symmetric four-quark interaction (nonlocal or local NJL type) 
the Goldstone particles and other mesons appear as the poles 
in the quark-antiquark scattering matrix due to the summa¬ 
tion of infinite number of diagrams [20, 27, 32, 37, 38, 40]. 
In these diagrams, the quark and antiquark interact via the 
four-quark interaction. On the other hand, in the box dia¬ 
gram (Fig. 2), the quark and antiquark do not interact be¬ 
tween each other and thus it is separated from the set of dia¬ 
grams producing mesons as bound states. It means in partic¬ 
ular, that in these approaches there are no double-counting 
effects. On the other hand in the framework NJL model it 
was shown that these two type of contributions, e.g. box ans 
bound state one, are necessary for the correct description of 
such processes as pion polarizability [52] or /r/r-scattering 
[53] and omitting one of these contribution will lead to large 
breaking of chiral symmetry. 

To above we can add that from the quark-hadron duality 
arguments, the quark loop (as for the two-point correlator as 
well for the four-point correlator) represents the contribution 
of the continuum of excited hadronic states. In the language 
of the spectral densities, the model calculations correspond 
to the model of the spectral density saturated by the lowest 
hadronic resonance plus the excited hadronic state contin¬ 
uum. The first part is for the meson-exchange diagrams, and 
the latter for the quark loop. It is the quark loop (continuum) 
provides the correct large photon momentum QCD asymp¬ 
totics for the Adler function, three- and four- point correla¬ 
tors. 

With the Feynman rules for the dynamical quark propa¬ 
gator (24) and the quark-photon vertices (29), (A.2), (A. 3), 
and (A.4), the gauge invariant set of diagrams describing the 
polarization tensor (< 7 i +<l 2 ),k+qi, —k) due to 

the dynamical quark loop contribution is given in Fig. 2. 

4 The results 

For the numerical estimates, the SU (2)- and SU (3)- versions 
of the N^QM model are used. In order to check the model 
dependence of the final results, we also perform calculations 
for different sets of model parameters. 



1.5 0.0 


Fig. 3 The 3D density p(Qi, Q 2 ) defined in Eqs. (20,21). 

In the SU (2) model, the same scheme of fixing the model 
parameters as in [38, 39] is applied: fitting the parameters A 
and m c by the physical values of the : r° mass and the K° —> 
yy decay width, and varying mo in the region 150 — 400 
MeV. For estimation of fl p LbL and its error, we use the re¬ 
gion for mo from 200 to 350 MeV. 

For the SU (3) version of the model, it is necessary to fix 
two more parameters: the current and dynamical masses of 
the strange quark. We suggest to fix them by fitting the K° 
mass and obtaining more or less reasonable values for the 77 
meson mass and the T] —> yy decay width. The main problem 
here is that the lowest value for the nonstrange dynamical 
mass mo is 240 MeV, because at lower mo the 77 meson 
becomes unstable within the model approach. 

Additionally, in order to show that the different schemes 
of parameter fixing will lead to similar results for a^ LbL , 
we calculate this quantity for the model ( 22 ) with param¬ 
eters taken from [54] for the Gaussian (G/ — G/y) and the 
Lorentzian (L/ —Ljy) nonlocal form factors. The authors of 
[54] have used other scheme of parameter fixing. Namely, 
the value of light current quark mass is fixed (8.5 MeV for 
G/ - Gnu 7.5 MeV for G IV , 4.0 MeV for L, - G m , and 3.5 
MeV for L/y). The other parameters are fitted in order to 
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Q^GeV] 

Fig. 4 The 2D slice of the density p{Q\,Qi) at Q 2 = Q\. Different 
curves correspond to the contributions of topologically different sets 
of diagrams drawn in Fig. 2. The contribution of the box diagram with 
the local vertices, Fig. 2a, is the dot (olive) line(Loc); the box diagram, 
Fig. 2a, with the nonlocal parts of the vertices is the dash (red) line 
(NLi); the triangle, Fig. 2b, and loop, Fig. 2c, diagrams with the two- 
photon vertices is the dash-dot (blue) line (NLi); the loop with the 
three-photon vertex, Fig. 2d, is the dot-dot (magenta) line (NL 3 ); the 
loop with the four-photon vertex, Fig. 2e, is the dash-dot-dot (green) 
line (NL 4 ); the sum of all contributions (Total) is the solid (black) line. 
At zero all contributions are finite. 

reproduce the values of the pion and kaon masses, the pion 
decay constant f n , and, alternatively, the 77 ' mass for sets G/, 
Giv, Lj, Ljv or the rj' —> yy decay width for sets G//, G///, 
Ln, Lijj. 

The important result, independent of the parameteriza- 
tions, is the behavior of the density p HLhL (Qi,Q 2 ), shown in 
Fig. 3. One can see, that p HLbL (Qi,Q 2 ) is zero at the edges 
(6i=0 or Qo = 0) and is concentrated in the low-energy 
region 7 (<2i ~ Q 2 ~ 300 MeV) providing the dominant con¬ 
tribution to a[, ILbL . This behavior at the edges appears to be 
due to cancelations of contributions from different diagrams 
of Fig. 2. 

In Fig. 4, the slice of p HLbL (Qi,Q 2 ) in the diagonal di¬ 
rection Q 2 = Q\ is presented together with the partial con¬ 
tributions from the diagrams of different topology. One can 
see, that the p HLbL ( 0 , 0 ) = 0 is due to a nontrivial cance¬ 
lation of different diagrams of Fig. 2. This important result 
is a consequence of gauge invariance and the spontaneous 
violation of the chiral symmetry, and represents the low en¬ 
ergy theorem analogous to the theorem for the Adler func¬ 
tion at zero momentum. Another interesting feature is, that 
the large Q\, Q 2 behavior is dominated by the box diagram 
with local vertices and quark propagators with momentum- 
independent masses in accordance with perturbative theory. 
All this is very important characteristics of the NjQM, in¬ 
terpolating the well-known results of the chiral perturbative 
theory at low momenta and the operator product expansion 

7 One should point out that the density for the mesonic exchanges has 

similar behavior. 


at large momenta. Earlier, similar results were obtained for 
the two-point [33, 34] and three-point [35] correlators. 

The numerical results for the value of flJJ LbL are given 
in the table and presented in Fig. 5 for the SU (2) and SU (3) 
models together with the result of C^QM [32] and DSE [40] 
calculations. The estimates for the partial contributions to 
a ttL bL 10 -iO) are ^ tj .0 contribution 5.01(0.37) [38], 
the sum of the contributions from 7 r°, 77 and 77 ' mesons 
5.85(0.87) [38], the scalar a, ao(980) and /o(980) mesons 
contribution 0.34(0.48) [3, 39], and the quark loop contri¬ 
bution is 11.0(0.9) [3]. In all cases we estimate the absolute 
value of the result and its error by calculating a |) LbL N zQ M 
for the space of model parameters fixed by above mentioned 
observables, except one, varying mo- Because in all cases 
the resulting curves (Fig. 5) are quite smooth, it gives to us a 
credit to point out rather small model errors (< 10 %) for the 
intermediate and final results. Thus our claim is that the total 
contribution obtained in the leading order in the 1 JN C expan¬ 
sion within the nonlocal chiral quark model is (see also [3]) 

rtpLbL’N^QM _ 16 . 8 ( 1-25 ). 10 — 10 (31) 

This value accounts for the spread of the results depending 
on reasonable variation of the model parameters and sensi¬ 
tivity to the different choice of the nonlocality shapes. Note, 
that as it was emphasized in [32], the results of these kind 
of calculations do not include the ’’systematic error” of the 
models. 

Comparing with other model calculations, we conclude 
that our results are quite close to the recent results obtained 
in [32, 40]. 8 It is no accidental. The most close to our model 
is the Dyson-Schwinger model used in [40]. The specific 
feature of both models is that the kernel of nonlocal inter¬ 
action is motivated by QCD. In [40] the kernel of the inter¬ 
action is generated by the nonperturbative gluon exchanges. 
In the N^QM the form of the kernel is motivated by the in- 
stanton vacuum models. The other difference between the 
NjQM and [40] is that, in a sense, the NjQM has a mini¬ 
mal structure (with respect to number of Lorentz structures 
for vertices, etc). Nevertheless, the predictions of the N)y QM 
for the different contributions to the muon g — 2 are in agree¬ 
ment with [40] within 10%. 

The constituent chiral quark model used in [32] corre¬ 
sponds to the local limit of the N^QM. This limit is achieved 
when the nonlocality parameter A goes to infinity, that means 
that the nonlocal form factors become constants: F(k 2 ,p 2 ) —> 
1. Taking this limit the N^QM becomes one-parametric one 
(only M q ) and we reproduce the M q dependence of quark 
box contribution to aJJ LbL shown in Fig. 13 of [32], What is 
more interesting and important is that, the M q dependencies 

s In earlier works [20, 27], the quark loop contribution was found one 
order less than in more recent calculations. To our opinion, one of the 
reason for that, is that in those models the photon-quark coupling is 
suppressed by the VMD form factors. 
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m d , u [MeV] m^pVleV] 


Fig. 5 (Left) The results for a{J LbL in the 5(7(2) model: the red dashed line is the total result, the green dotted line is the quark loop contribution 
and the magenta dash-dot-dot line is the K + a contribution. Thin vertical line indicates the region for estimation of aJ] LbL error band. (Right) The 
results for the black solid line is the 5(7(3)-result, the red dash line corresponds to the 5(7(2)-result, the blue dash-dotted line is the C%QM 

result [32], hatched region correspond to DSE result [40], 


of total contribution to aJJ LbL in [32] (Fig. 14) and in the 
N^QM have the same qualitative behavior and very close 
(with less than 10%) qualitatively. This is clear from Fig. 5. 

These facts are very pleasant for the phenomenology of 
the HLbL contributions to the muon, because it means that 
even starting from the models that differ in many details, the 
predictions are still very stable numerically. 


5 Conclusions 

In this paper, we have presented the results for the con¬ 
tribution of the dynamical quark loop mechanism for the 
light-by-light scattering to the muon anomalous magnetic 
moment within the nonlocal chiral quark model. In previ¬ 
ous works [3, 38, 39], we calculated the corresponding con¬ 
tributions due to the exchange by pseudoscalar and scalar 
mesons. The basis of our model calculations is the sponta¬ 
neous violation of the chiral symmetry in the model with 
the nonlocal four-fermion interaction and abelian gauge in¬ 
variance. The first leads to the generation of the momentum- 
dependent dynamical quark mass, and the latter ensures the 
fulfillment of the Ward-Takahashi identities with respect to 
the quark-photon interaction. 

In the present work, we derived the general expression 
for a bbL as the three-dimensional integral in modulus of the 
two photon momenta and the angle between them. The inte¬ 
gral is the convolution of the known kinematical factors and 
some projections of the four-photon polarization tensor. The 
latter is the subject of theoretical calculations. 

Since our model calculations of the hadronic contribu¬ 
tions are basically numerical, it is more convenient to present 
our results in terms of the density function p HLbL (Qi,Q 2 )- 


We observe some properties of this function that have model- 
independent character. Firstly, at zero momenta one has 
p hi, hi.('o o) — o in spite of the fact that the partial contri¬ 
butions of different diagrams are nonzero in this limit. This 
low-energy theorem is a direct consequence of the quark- 
photon gauge invariance and the spontaneous violation of 
the chiral symmetry. Secondly, at high momenta the density 
is saturated by the contribution from the box diagram with 
the local quark-photon vertices and local quark propagators 
in accordance with the perturbative theory. This is a conse¬ 
quence of the fact, that at small distances all nonperturbative 
nonlocal effects are washed out. Thirdly, with the model pa¬ 
rameters chosen, the p IILbL ((7|. Qo ) is concentrated in the 
region gi ~ ^ 300 MeV, which is a typical scale for 

light hadrons. 

Summarizing the results of the present and previous works 
[3, 38, 39], we get the total hadronic contribution to aJJ LbL 
within the N^QM in the leading order in the 1 /N c expan¬ 
sion. The total result is given in Eq. (31). To estimate the 
uncertainty of this result, we vary some of the model pa¬ 
rameters in physically reasonable interval and also study the 
sensitivity of the result with respect to different model pa- 
rameterizations. In this sense, the error in Eq. (31) is a con¬ 
servative one. 

If we add the result (31) to all other known contributions 
of the standard model to a u , (3)-(7), we get that the differ¬ 
ence between experiment (2) and theory is 
fl BNL,CODATA _ a SM = jg 73 x IQ- 10 , (32 ) 

which corresponds to 2.43(7. If one uses the hadronic vac¬ 
uum polarization contribution from the T hadronic decays 
instead of e + e~ data 

^hvplo-t = 701 5(4 7) x 10 -io [15]; (33) 


















the difference decreases to 18.44 x 10 10 (2.23(7) for the 
case of aJJ LbL from (8) [16] and to 12.14 x 10~ 10 (1.53(7) in 
our model (31). 

Clearly, a further reduction of both the experimental and 
theoretical uncertainties is necessary. On the theoretical side, 
the calculation of the still badly known hadronic light-by- 
light contributions in the next-to-leading order in the 1 /N c 
expansion (the pion and kaon loops) and extension of the 
model by including heavier vector and axial-vector mesons 
is the next goal. The contribution of these effects and the 
model error induced by them are not included in the result 
(31). Preliminary studies [20, 22] show that these contribu¬ 
tions are one order smaller than the pseudoscalar exchanges 
and the quark loop contributions. However, the interesting 
point that inclusion of vector channel can strongly suppress 
contribution from the quark loop due photon-vector meson 
exchange which lead to appearance in each photon vertex 
additional VMD-like factor. This was found in local NJL 
model [27] and should be carefully investigated in the non¬ 
local one. 

Work in this direction is now in progress, and we hope 
to report its results in the near future. 
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Appendix A: Nonlocal multi-photon vertices 

Let us introduce the finite-difference derivatives 


At)/ M _ f{a + b)-f(b) 


/ 1 ( a >b) = 


(i a + b)~ — b 2 


(A.l) 


f in+1 H a Abi},b h b 2 ) = 


f (n) (a, {£>/},£>i) — (a,{bj},b 2 ) 
(a + bi) 2 - (a + b 2 ) 2 


where n = 1,2,... . Then, the quark-antiquark vertex with 
the two-photon insertions (Fig. lb) is [50] 


rffi (quqi) = 2gp V »t (1) (k,k r ) + 

(k + ki (*t + k') y mP) (*, h , H) + (A.2) 

(k + k 2 ) v (k 2 + k l )^m {2] (k,k 2 ,k r ) . 

Here and below, k is the momentum of the incoming quark, 
k' is the momentum of the outgoing quark, q-, are the mo¬ 
menta of the incoming photons, and k\ = k + qi, kjj j, = 
k -\-qi ~\~qj + + qk- 


The quark-three-photon vertex (Fig. lc) is [51] 


J]ivp(<7i,<72,<73) = - [2g ll v{kn + k , )pm (2 \kM 2 ^) 
+ 2gp V (k + k 3 ) p iii (2 \k } k 3} k') 


+ (k + ki)n(ki +k n ) v (kx2+k')piiP\k 1 k Xl ki 2 ,k') 

+ (k + ki)p(ki 3 +l^) v (ki +k\ 2 )pinP'\k,k\,k\ 2 ,k')\ 

+ [1^3,p^p] + [2^3,v^p], (A.3) 


The quark-four-photon vertex (Fig. Id) takes the form 


rW 

1 l~LVp Z 


( 91 , <72, <73, < 74 ) 


+ 4gpvgTpm {2 ' , (k 1 k l2 ,k') 


+ 4gp V g T pm [2 \k,k 24 ,k') 


+ 2gp V (^(k+k 2 )p(k 3 , + k 3 , 4 ) t m {i \k 1 k i ,k 2 , 4 ,k') 

+ {k + k 2 ) p (k m + k') x m (3) (k, k 3 ,ki 2 i,k') 

+ (kn + km)p{km + k , ) T m( 3 ') (k,kn,km,k') 

+ {ki 24 + k')p {k\ 2 + &i24)iV 3) (k, kn,ku 4 ,k') 

+ {k\24 + k') p (k + k 4 ) x m^'> (k, k 4 ,ki24,k') 

+ (k4 + k 2 4) p (k + k4) T mP\k,k4,k 2 ,4,k!) s j 

+ 2g zp ({k+ki)p(k\ +k l2 ) v m {i \k 1 k u ki 2 ,k') 

+ (k + k 2 ) v (k 2 +ki2}pm^\k,k 2 ,ki2,k!) 

+ (£34 + ^234 ) V (&234 + &1 234 ) p m ^ * (k, £ 34 , ^234, k') 
+ (£34 + ki 2 4 )p (&134 + ^ 1234 ) yirP\k, &34, k\ 3 4 ,k!) 
+ (k + kx)p(ki34 + ki234)v>n^ {k,kx,ki34,k') 

+ (k + k 2 ) v (k 22 4 + fci234)/pra (3) (k, k 2 , k 23 4 , k')^j 


(A.4) 


+ (k + ki)p(k\ +k\ 2 ) v [k\ 2 + ki 22 )p(ki 22 + k') x x 

x m^\k,ki,ki2,ki23,k') 

+ (k + ki)p(ki +k\ 2 ) v (ki 2 4 + k\ 2 34 )p{k \2 + k\ 2 4 ) x x 

x (k,k\,ki2,k\2A,k!) 

+ (k 2 + k\ 2 )p (k + k 2 ) v (k\ 2 + ki 22 )p(ki 22 + £1234)1 x 

x m ^ {k,k 2 ,ki 2 ,ki 22 ,k!) 

+ (k 2 + k\ 2 )p(k + k 2 ) v (k\ 2 4 + k 123 4)p(k l2 + k 124 ) x x 
x 111 {k,k 2 ,ki 2 ,ki 24 ,kf) 

+ (k 23 + ki 23 )p{k +k 2 )v{k 2 + k 23 )p (£123 + £1234)1 x 
x m {k,k 2 ,k 23 ,k\ 23 ,k!) 

+ [k 2 4 + k\ 2 4)p{k + k 2 ) v {k\ 2 4 + k\ 23 4) p(k 2 + k 2 4) x x 
x m^\k,k2MAMu,k!) 

+ (k 23 4 J rk\ 23 4 )p(k + k 2 )v(k 2 4 + k 23 4 )p{k 2 +k 2 4 ) x x 

x (k,k 2 ,k 2 4 ,k 234 ,k?) 

+ (£234 + £1234) p (£ + £2) v (£2 + £23 )p (£23 + ^234) t x 


x m(k, k 2 , k 23 , k 23 4 , k') 
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+ [2^4,v^T] + [2^3,v^p]. 
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model parameters 



fl bbL in 

io- 10 



No. 

m d ,u 

m c ,u 

m d ,s 

m c ,s 

A 

total 

total 

u + d 

4’ quark 

w o,„ 

n + ??'+ 


MeV 

MeV 

MeV 

MeV 

MeV 

517(3) 

517(2) 

quark loop 

loop 


ao + fo 

1 

150 

0.33 



6786.6 


21.74 

16.95 


4.79 


2 

160 

0.59 



4890.7 


20.49 

15.38 


5.11 


3 

170 

0.93 



3768.9 


19.45 

14.10 


5.35 


4 

180 

1.33 



3049.2 


18.60 

13.10 


5.50 


5 

190 

1.78 



2557.6 


17.89 

12.29 


5.60 


6 

200 

2.27 



2204.9 


17.30 

11.67 


5.64 


7 

210 

2.79 



1941.5 


16.82 

11.18 


5.64 


8 

220 

3.34 



1738.2 


16.42 

10.80 


5.62 


9 

230 

3.90 



1577.1 


16.09 

10.52 


5.58 


10 

240 

4.47 



1446.3 


15.83 

10.31 


5.52 


11 

240 

4.47 

339.5 

133.7 

1446.3 

16.53 

15.89 

10.31 

0.234 

5.58 

0.404 

12 

250 

5.06 



1338.2 


15.61 

10.16 


5.46 


13 

250 

5.06 

347.1 

148.4 

1338.2 

16.21 

15.64 

10.16 

0.229 

5.49 

0.337 

14 

260 

5.65 



1247.2 


15.45 

10.05 


5.39 


15 

270 

6.25 



1169.6 


15.31 

9.99 


5.33 


16 

280 

6.86 



1102.5 


15.23 

9.97 


5.26 


17 

280 

6.86 

387.4 

193.4 

1102.5 

15.91 

15.22 

9.97 

0.209 

5.25 

0.480 

18 

290 

7.48 



1043.9 


15.14 

9.94 


5.19 


19 

300 

8.09 



992.2 


15.08 

9.96 


5.12 


20 

305 

8.41 

413.7 

231.9 

968.6 

15.78 

15.05 

9.97 

0.193 

5.08 

0.543 

21 

310 

8.72 



946.2 


15.04 

9.98 


5.06 


22 

320 

9.35 



904.9 


15.02 

10.02 


5.00 


23 

320 

9.35 

428 

255.7 

904.9 

15.78 

15.02 

10.03 

0.182 

4.99 

0.577 

24 

330 

9.99 



867.7 


15.01 

10.08 


4.94 


25 

340 

10.63 



833.8 


15.02 

10.14 


4.88 


26 

350 

11.29 



802.8 


15.03 

10.21 


4.82 


27 

350 

11.29 

451.2 

305.9 

802.8 

15.80 

15.05 

10.21 

0.159 

4.84 

0.592 

28 

360 

11.95 



774.4 


15.07 

10.30 


4.77 


29 

370 

12.62 



748.1 


15.10 

10.38 


4.72 


30 

380 

13.30 



723.8 


15.14 

10.47 


4.67 


31 

390 

13.99 



701.1 


15.19 

10.56 


4.63 


32 

400 

14.69 



679.8 


15.24 

10.65 


4.58 


G, 

304.5 

8.50 

All 

223 

1002.7 

15.67 

14.67 

9.71 

0.192 

4.95 

0.810 

G // 

304.5 

8.50 

439 

223 

1002.7 

15.93 

14.67 

9.71 

0.190 

4.95 

1.070 

Gin 

304.5 

8.50 

422 

223 

1002.7 

15.57 

14.67 

9.71 

0.193 

4.95 

0.707 

G iv 

287.5 

7.50 

408 

199 

1086.1 

15.75 

14.81 

9.75 

0.202 

5.06 

0.738 

L, 

295 

4.00 

450 

112 

1013 

15.58 

14.84 

9.61 

0.241 

5.23 

0.503 

Lji 

295 

4.00 

505 

110 

1013 

16.37 

14.84 

9.61 

0.222 

5.23 

1.311 

L iii 

296 

4.00 

457 

112 

1013 

15.61 

14.80 

9.57 

0.238 

5.22 

0.578 

L iv 

277.5 

3.50 

418 

100 

1110 

15.76 

15.08 

9.74 

0.251 

5.33 

0.435 


Table 1 The contribution to the muon AMM a bbL for different sets of model parameters. The model parameters G/_/y and L j_jy are taken from 
[54], The difference between a bbL for the set of model parameters with same u quark mass in SU (2) and SU (3) models is due to mixing of <7 meson 
with /o meson is SU (3) case. In order to extrapolate the SU (3)-result to lower quark masses, we find that the maximal value of the difference 
between the total SU (2)- and SU (3)-results is 0.77. We add this number to the value of the SU (2)-result at 200 MeV. 















